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APPLICATIONS OF THE CONTRAVARIANT FORM OF THE NAVIER-STOKES EQUATIONS 


Theodore Katsanis 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

The contravariant Navier-Stokes equations in weak conservation form are 
well suited to certain fluid flow analysis problems. Three-dimensional con- 
travariant momentum equations may be used to obtain Navier-Stokes equations in 
weak conservation form on a nonplanar two-dimensional surface with varying 
streamsheet thickness. Thus a three-dimensional flow can be simulated with 
two-dimensional equations to obtain a quasi-three-dimensional solution for 
viscous flow. When the Navier-Stokes equations on the two-dimensional non- 
planar surface are transformed to a generalized body-fitted mesh coordinate 
system, the resulting equations are similar to the equations for a body-fitted 
mesh coordinate system on the Euclidean plane. 

Contravariant momentum components are also useful for analyzing compres- 
sible, three-dimensional viscous flow through an internal duct by parabolic 
marching. This type of flow can be efficiently analyzed by parabolic marching 
methods, where the streamwise momentum equation is uncoupled from the two 
crossflow momentum equations. This can be done, even for ducts with a large 
amount of turning, if the Navier-Stokes equations are written with contravar- 
iant components. 


INTRODUCTION 

The Navier-Stokes equations are often transformed to special coordinate 
systems for computational purposes. However, if it is desired to have the 
equations in conservation form, it is generally necessary to use Cartesian mo- 
mentum components. On the other hand, there are cases where it is useful to 
use the contravariant momentum components in weak conservation form. Two ex- 
amples are considered here: one example is a nonplanar flow surface embedded 

in a three-dimensional flow field; the other example is three-dimensional flow 
in a turning duct. Both examples have applications in turbomachinery. 

There are many instances of two-dimensional flow where the flow is not 
planar; i.e., either the two-dimensional surface is curved or the streamsheet 
thickness is not constant (e.g., axisymmetric flow). In such cases the three- 
dimensional Navier-Stokes equations may be reduced to two dimensions. A gen- 
eral method for doing this is to transform the three-dimensional equation to a 
new coordinate system such that two coordinate lines lie on the surface of in- 
terest, with the third coordinate normal to the surface. The momentum equa- 
tions are then combined so as to obtain contravariant momentum components. 

Thus two components of momentum lie on the surface and the third component is 
omitted; the equation will involve only the two contravariant components. As 
an illustration, equations are presented here for viscous flow on a rotating 
surface of revolution. A flow solution on such a surface would be useful for 
analyzing flow through a turbomachine. 
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Ducts of complex geometry may be analyzed efficiently if the flow is of a 
"boundary layer" type; i.e., flow with a predominant direction (no reverse 
flow), and negligible streamwise diffusion. This kind of flow can be analyzed 
by parabolic marching methods. Patankar and Spalding (ref. 4), Briley 
(ref. 8), Roberts and Forester (ref. 7), Moore and Moore (ref. 6), and others 
have calculated three-dimensional viscous flow in ducts by means of parabolic 
space marching. For this purpose the momentum equations are uncoupled. The 
streamwise momentum equation is used to advance streamwise velocity components, 
and cross-stream momentum equations are used to advance the cross-stream veloc- 
ity components. If a body-fitted coordinate system is used, the contravariant 
components of the momentum equation (with contravariant velocity components as 
dependent variables) provide a natural method of uncoupling the streamwise mo- 
mentum equation from the cross-stream momentum equation. As an example, equa- 
tions are presented for flow through a rotating passage, e.g., a blade passage 
through a turbomachine blade row. 


NAVIER-STOKES EQUATIONS 

The two-dimensional, planar, nonsteady Navier-Stokes equations in conser- 
vation form are 
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Here p is density, V is velocity, p is pressure, e is total internal en- 
ergy per unit volume, T is temperature, and t is shear stress tensor. The 
x and y subscripts refer to the x and y coordinates. For numerical so- 
lutions, the equation may be transformed to an arbitrary ^ - n coordinate 
system (ref. 1). These are usually body-fitted coordinates, such that constant 
values of c or n correspond to the body surface. The transformed equations 
in conservation form are 
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The superscripts 5 and n refer to contravariant velocity components, 
and J is the Jacobian of the transformation. The thin-layer approximation 
is used here (viscous derivatives in the streamwise direction are neglected). 
Note that the momentum components (the second and third equations) are still x 
and y components. If £ and n (either contravariant or covariant) momen- 
tum components are used the equation cannot be written in conservation form if 
there is any curvature of the t; or n coordinate lines. 

It is desired to obtain two-dimensional Navier-Stokes equations for non- 
pi anar surfaces (i.e., quasi-three-dimensional). The desired equations can be 
obtained from the three-dimensional version of Eq. (2), written with contra- 
variant momentum components in a 5 , n, s coordinate system. The 5 and n 
coordinates lie in the desired nonplanar surface, and 5 is chosen in a spe- 
cial way to be orthogonal to the surface. It is assumed that the chosen non- 
planar surface is a stream surface, i.e., the contravariant velocity component 
V? =0. When this is done the third momentum equation can be eliminated and 
the two-dimensional Navier-Stokes equation for flow on a two-dimensional sur- 
face is obtained. Rather than give the general equation, an example will be 
given which is of interest for turbomachinery flow analysis. 


NAVIER-STOKES EQUATIONS ON A ROTATING SURFACE OF REVOLUTION 

For analysis of flow through turbomachinery blade rows it is useful to 
consider flow on a blade-to-blade surface as shown in figure 1. The coordi- 
nates are m (physical distance) and e (radians). In addition the stream- 
sheet thickness h must be specified, as shown in figure 2. The thickness 
h and the radius r are both functions of m only. The three-dimensional 
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contravariant momentum equations are obtained in m, e, n coordinates, where 
n is normal to the blade-to-blade surface. The n coordinate is defined such 
that the metric coefficient is 1/h, so that not only is the normal velocity 
zero, but also the derivative in the n direction is zero. The surface of 
revolution is rotating at a constant rotational velocity u . The relative 
velocity is W; m and e subscripts will denote physical (not covariant) 
velocities. 


The Navier-Stokes equation obtained is 
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= 1/rh is the Jacobian and m and e momentum components are used. 
Except for the source term, K', Eq. (3) is similar to Eq. (1). Eq. (3) can now 

be transformed to body-fitted coordinates, similar to Eq. (2). The transformed 
equation is 

3 t q' + 3„E' + 3 F' = a S' + K* ( 4 ) 
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J = J1J2 is the Jacobian (x, y, z) + (5, n, n), J2 is the Jacobian 
( m » 0 ) * U , n), and m and e momentum components are used. Equation (4) 
uses m and e momentum components, and is similar to Eq. ( 2 ), except for 
the source term K‘ which is required because of the curvature of the m and 
0 coordinates. The techniques and computer codes for solving Eq. ( 2 ) (ref. 1 ) 
should be able to be extended to Eq. ( 4 ). 

THREE-DIMENSIONAL NAVIER-STOKES EQUATIONS FOR A ROTATING COORDINATE SYSTEM 

An example of a viscous, three-dimensional internal flow that can be ana- 
lyzed by parabolic space marching is flow through a turbomachinery blade row. 
The desired equations are obtained by transforming the Navier— Stokes equations 
in cylindrical (r, e, z) coordinates to a body-fitted (e, n , c) coordinate sys- 
tem. Then the contravariant components of the momentum equations are obtained. 
These equations involve only the contravariant velocity components. (The co- 
variant components of momentum have a simpler form, but involve both covariant 
and contravariant velocity components (ref. 5 ).) This permits the streamwise 
momentum calculation to be uncoupled from the cross-stream momentum calcula- 
tion. Continuity is satisfied indirectly by modifying an inviscid pressure 
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field to correct for viscous effects. Modifying the streamwise pressure gra- 
dient will affect streamwise velocity components so that global continuity is 
satisfied. Modifying the pressure over a cross section will affect the cross- 
stream velocities so that local continuity can be satisfied. The steady-state 
three-dimensional Navier-Stokes equations for a rotating coordinate system 
(transformed to a body-fitted (5, n, c) coordinate system, and with contravar- 
iant momentum components) are: 
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The flow is in the 5 direction, and n, c are the cross-stream directions. 
The thin layer approximation is used. 


CONCLUDING REMARKS 

Two examples of applications of the contravariant Navier-Stokes equations 
have been presented. 

The first example shows that the contravariant Navier-Stokes equations 
can be used to obtain two-dimensional equations for a rotating surface of 
revolution with varying streamsheet thickness. The equations are similar in 
form to the conservation-law form of the Cartesian two-dimensional Navier- 
Stokes equations. The methods already developed for solving the Cartesian 
two-dimensional equations can thus be used for analyzing flow on a blade-to- 
blade surface of a turbomachine. Other applications for nonplanar surfaces 
seem possible. 

The second example is for three-dimensional viscous flow through a rota- 
ting duct. The contravariant form of the Navier-Stokes equations provides a 
method for uncoupling the streamwise momentum from the cross-stream momentum. 
This makes it easier to use parabolic marching methods for turbomachinery 
blade passages with a large amount of turning. A similar technique could be 
used with other types of internal passages. 
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